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Abstract 

Topological simplification of scalar and vector fields is well- 
established as an effective method for analysing and visualising 
complex data sets. For multi-field data, topological analysis re¬ 
quires simultaneous advances both mathematically and computa¬ 
tionally. We propose a robust multivariate topology simplifica¬ 
tion method based on “lip”-pruning from the Reeb Space. Math¬ 
ematically, we show that the projection of the Jacobi Set of mul¬ 
tivariate data into the Reeb Space produces a Jacobi Structure 
that separates the Reeb Space into simple components. We also 
show that the dual graph of these components gives rise to a Reeb 
Skeleton that has properties similar to the scalar contour tree and 
Reeb Graph, for topologically simple domains. We then intro¬ 
duce a range measure to give a scaling-invariant total ordering 
of the components or features that can be used for simplifica¬ 
tion. Computationally, we show how to compute Jacobi Struc¬ 
ture, Reeb Skeleton, Range and Geometric Measures in the Joint 
Contour Net (an approximation of the Reeb Space) and that these 
can be used for visualisation similar to the contour tree or Reeb 
Graph. 

keyword Simplification, Topology, Multi-Field, Reeb Space, 
Joint Contour Net, Multi-Dimensional Reeb Graph, Reeb Skele¬ 
ton 

1 Introduction 

Scientific data is often complex in nature and difficult to vi¬ 
sualise. As a result, analytic tools have become increasingly 
prominent in scientific visualisation, and in particular topolog¬ 
ical analysis. While earlier work dealt primarily with scalar 
data (331 lam, multivariate topological analysis in the form of 
the Reeb Space HSl [T9ll has started to become feasible using a 
quantised approximation called the Joint Contour Net (JCN) Q . 

Prior experience in scalar and vector topology shows that sim¬ 
plification of topological structures is required, as real data sets 
are often noisy and complex. Although most of the work required 
is practical and algorithmic in nature, mathematical formalisms 
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are also needed, in this case based on fiber analysis, in the same 
way that Reeb Graphs and contour trees rely on Morse theory. 
This paper therefore: 

1. Clarifies relationships between the Reeb Space of a multi¬ 
variate map /, the Jacobi Set of /, and fiber topology, 

2. Introduces the Jacobi Structure in the Reeb Space that de¬ 
composes the Reeb Space into regular and singular compo¬ 
nents equivalent to edges and vertices in the Reeb Graph, 
then reduces it further to a Reeb Skeleton, 

3. Proves that Reeb Spaces for topologically simple domains 
have simple structures with properties analogous to proper¬ 
ties of the contour tree, allowing lip-pruning based simplifi¬ 
cation, 

4. Introduces the range measure and other geometric mea¬ 
sures for a total ordering of regular components of the Reeb 
Space, 

5. Describes an algorithm that extracts the Jacobi Structure 
from the Joint Contour Net using a Multi-dimensional Reeb 
Graph (MDRG) and computes the Reeb Skeleton, and 

6. Simplifies the Reeb Skeleton and the corresponding Reeb 
Space computing the range and other geometric measures 
using the Joint Contour Net. 

To clarify the relationships between the newly introduced data- 
structures in the current paper, note that the JCN is an approxi¬ 
mation of the Reeb Space. We compute a MDRG from the JCN. 
The critical nodes of the MDRG form the Jacobi Structure of the 
JCN. The Jacobi Structure then separates the JCN into regular 
and singular components. The dual graph of such components 
gives a Reeb Skeleton which is used in the multivariate topology 
simplification. 

As a result, much of this paper addresses the theoretical ma¬ 
chinery for simplification of the Reeb Space and its approxima¬ 
tion, the Joint Contour Net. Section [2| reviews relevant back¬ 
ground material on simplification, followed by a more detailed 
review of the fiber topology, Jacobi Set and Reeb Space in Sec¬ 
tion [S] Section [4] provides theoretical analysis and results needed 
for the lip-simplification of the Reeb Space. For simple domains, 
the Reeb Space can have detachable (lip) components: this is 
used in Section [5] to generalise leaf-pruning simplification from 
the contour tree to the Reeb Space. Once this has been done, we 
introduce a range persistence and other geometric measures to 
govern the simplification process. 

In Section [6l we give an algorithm for simplifying the Joint 
Contour Net (an approximation of the Reeb Space). We start 
by building a hierarchical structure called the Multi-Dimensional 
Reeb Graph (MDRG) that captures the Jacobi Structure of the 
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Joint Contour Net, and then show how to reduce the JCN to a 
Reeb Skeleton - a graph with properties similar to a contour tree. 
In Section[71 we illustrate these reductions first with analytic data 
where the correct solution is known a priori, then for a real data 
from the nuclear physics. As part of this, we provide performance 
figures and other implementation details in Section [71 then draw 
conclusions and lay out a road map for further work in Section[8l 

2 Previous Work 

Topology-based simplification aims to reduce the topological 
complexity of the underlying data. There are different ways to 
measure such topological complexity depending on the nature 
of the underlying data. Here we mention some well-known ap¬ 
proaches from the literature for measuring the topological com¬ 
plexity and their simplification procedure. 

Scalar Field Simplification. 

The topological complexity of the scalar field data is measured 
in terms of the number of critical points and their connectivities - 
captured by its Reeb Graph or contour tree. Another way to cap¬ 
ture the topological complexity of the scalar field is by computing 
the Morse-Smale complex of the corresponding gradient field. 
Therefore, the topological simplification in this case is driven by 
reducing the number of critical points via simplification of the 
Reeb Graph/ contour tree or the Morse-Smale complex. Carr et 
al. 0 describe a method for associating local geometric mea¬ 
sures such as the surface area and the contained volume of con¬ 
tours with the contour tree and then simplifying the contour tree 
by suppressing the minor topological features of the data. Note 
that a feature is any prominent or distinctive part or quality that 
characterises the data and topological features captures the topo¬ 
logical phenomena of the underlying data. Wood et al. 1(5^ give 
a Reeb Graph based simplification strategy for removing the ex¬ 
cess topology created by unwanted handles in an isosurface using 
a measure for computing the handle-size in the isosurface and as¬ 
sociating them with the loops of the Reeb Graph. Gyulassy et 
al. |[24l describe a technique for simplifying a three-dimensional 
scalar field by repeatedly removing pair of critical points from the 
Morse-Smale complex of its gradient field, by repeated applica¬ 
tion of a critical-point simplification operation. Mathematically, 
the simplification of “lips” proposed in this paper is a direct gen¬ 
eralization of this idea (for scalar fields) to multi-fields. Luo et al 
|[37]| describe a method for computing and simplifying gradients 
and critical points of a function from a point cloud. Tierny et al. 
||5^ present a combinatorial algorithm for simplifying the topol¬ 
ogy of a scalar field on a surface by approximating with a simpler 
scalar field having a subset of critical points of the given field, 
while guaranteeing a small error distance between the fields. 

The topological complexity of a point cloud data can be mea¬ 
sured by its homology. For a point cloud data in M? this is ex¬ 
pressed by the topological invariants, such as the Betti numbers 
corresponding to a simplicial complex of the point cloud - de¬ 
noted by Po (number of connected components). Pi (number of 
tunnels or 1-dimensional holes) and P 2 (number of voids or 2- 
dimensional holes). The i-th Betti number represents the rank of 
the i-th homology group (i = 0,1,2). Edelsbrunner et al. ll20t in¬ 


troduce the idea of persistence homology for the topological sim¬ 
plification of a point cloud by reducing the Betti numbers using 
a filtration technique. Cohen-Steiner et al. ITTll extend the per¬ 
sistence diagram for scalar functions on topological spaces and 
analyze its stability. 

Mesh Simplification. 

Mesh-simplification is well-known in the computational geome¬ 
try and graphics community. Topological complexity of a mesh 
can be determined by its genus. Guskov et al. 1231 remove the 
unnecessary topological noise from meshes of laser scanner data 
by reducing their genera. Nooruddin et al. Il40ll give a voxel- 
based simplification and repair method of polygonal models us¬ 
ing a volumetric morphological operation. Ni et al. |[39l gen¬ 
erate a fair Morse function for extracting the topological struc¬ 
ture of a surface mesh by user-controlled number and configu¬ 
ration of critical points. Hoppe et al. |[30l describe a energy- 
minimization technique for generating an optimal mesh by re¬ 
ducing the number of vertices from a given mesh. Also Hoppe et 
al. lEgm give a new progressive mesh representation, a new 
scheme for storing and transmitting arbitrary triangle meshes, 
and their simplification technique. Chiang et al. ll33l describe a 
technique of progressive simplification of tetrahedral meshes pre¬ 
serving isosurface topologies. Their method works in two stages 
- first they segment the volume data into topological-equivalence 
regions and in the second step they simplify each topological- 
equivalence region independently by edge collapsing, preserving 
the iso-surface topologies. There are many cost-driven methods 
of mesh-simplification (in the literature) which attempt to mea¬ 
sure only the cost of each individual edge collapse and the entire 
simplification process is considered as a sequence of steps of in¬ 
creasing cost (El [351 [33 ED. 

Vector Field Simplification. 

Topology based methods for vector field simplification are based 
on the idea of singularity pair cancellation to reduce the number 
of singularities and thus the topological complexity. This method 
iteratively eliminates suitable pairs of singularities with opposite 
Poincare-Hopf indices so that total sum of the indices remain in¬ 
variant to keep the global structure of the field the same. This idea 
has been exploited in ||56l[59l[45l. There are also non-topology 
based methods for vector-field simplification which are mainly 
based on smoothing operations. Smoothing operations reduce 
vector and tensor-field complexity and remove large percentage 
of singularities. Polthier et al. 1431 apply Laplacian smoothing 
on the potential of a vector-field. Tong et al. |[55l decompose a 
vector field into three components: curl free, divergence free and 
harmonic. Each component is smoothed individually and results 
are summed to obtain simplified vector field. 

Multi-Field Simplification. 

To the best of our knowledge, until now there is no prior work 
on topology-based simplification of general multi-field data. 
All those techniques, cited so far, for simplifying scalar fields, 
meshes and vector fields are not directly applicable in case of 
multi-fields, mainly because the computation of the equivalent 
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tools such as, Jacobi Set ca, Reeb Space M are not well- 
developed. A generalization of the persistence homology is 
proven to be difficult for the multi-fields ||4l . However, few at¬ 
tempts have been made for simplifying the Jacobi Sets in restric¬ 
tive cases. Snyder et al. im give two metrics for measuring per¬ 
sistence of the Jacobi Sets. Bremer et al. a describe a method 
for noise removal from the Jacobi Sets of time varying data. 
Suthambhara et al. |[52ll give a technique for the Jacobi Set sim¬ 
plification of bivariate fields based on simplification of the Reeb 
Graphs of their comparison measures. Huettenberger et al. pro¬ 
pose multi-field simplification method using Pareto sets ll3n[32ll . 
However, these methods lack mathematical justification for sim¬ 
plifying the corresponding input multi-fields and work mostly for 
bivariate data. In a similar context, Bhatia et al. O provide a sim¬ 
plification method by generalising the critical point cancellation 
of scalar functions to the Jacobi Sets in two dimensional domains. 

However, current research shows that the Jacobi Sets are un¬ 
able to capture the actual topological changes of multi-fields, 
instead one should consider their Reeb Spaces, introduced in 
M- Recently, Multi-Dimensional Reeb Graphs ifTOl and Lay¬ 
ered Reeb Graphs llSll have been introduced from two different 
perspectives to extend the Reeb Graph for multi-fields. In the cur¬ 
rent paper, we use the recently introduced Jacobi Structure m 
to separate the Reeb Space into regular and singular components. 
Thus we obtain a dual Reeb Skeleton corresponding to the Reeb 
Space. Our simplification strategy is based on simplifying this 
Reeb Skeleton by associating different measures with the nodes 
of the Reeb Skeleton. 

3 Necessary Background 

Over the last two decades, scalar topology has been used to 
support scientific data analysis and visualization, in particular 
through the use of the Reeb Graph and its specialisation, the con¬ 
tour tree lETiiaEalnisD. The subject of multi-field topology in 
data analysis is rather new. In this section we briefly describe the 
multi-field topological analysis and existing tools for capturing 
them, viz. the Jacobi Set and the Reeb Space. 


Table 1: Important Notations 


Notation 

Name 

W/ 

Reeb Space 

% 

Reeb Skeleton 

If 

Jacobi Set 

li 

Boundary Jacobi Set 

Ff 

Interior Jacobi Set 

3/ 

Jacobi Structure 

JCN{f,mQ) 

Joint Contour Net with quantization level mg 

% 

Reeb Graph 

My 

Multi-dimensional Reeb Graph 


Multi-Field Analysis. 

A multi-field on a J-manifold X (C with r component scalar 
fields /;•: X ^ M (/ = 1, ..., r) is a map / = (/i, / 2 , ..., /r) : X ^ 


W. Table [T] shows the notations used to denote various structures 
corresponding to a multi-field / in the current paper. 

In differential topology, / is considered to be a smooth map 
when all its partial derivatives of any order are continuous. A 
point X G X is called a singular point (or critical point) of / if 
the rank of its differential map dfx is strictly less than min{ J, r} 
where dfx is the r x J matrix whose rows are the gradients of fi to 
fr at X. And the corresponding value /(x) = c = (ci , C 2 , ..., 
in W is a singular value. Otherwise if the rank of the differential 
map dfx is min{ r} then x is called a regular point and a point 
y G is a regular value if f~^{y) does not contain a singular 
point. 

The inverse image of the map / corresponding to a value c G 
W, f~^{c) is called a fiber and each connected component of 
the fiber is called a fiber-component (491SH • In particular, for 
a scalar field these are known as the level set and the contour, 
respectively. The inverse image of a singular value is called a 
singular fiber and the inverse image of a regular value is called 
a regular fiber. If a fiber-component passes through a singular 
point, it is called a singular fiber-component. Otherwise, it is 
known as a regular fiber-component. Note that a singular fiber 
may contain a regular fiber-component. 

A continuous map is said to be proper if the pre-image of a 
compact set is always compact and it is said to be stable if its 
topological properties remain unchanged by small perturbations 
(281. Let / : X C ^ be a proper smooth map. Then, it 
is stable if and only if it satisfies the following local and global 
conditions. Around each singular point s, / is locally described as 
either (i) (t/, +y^): s is a definite fold point, or (ii) {u^x^ —y^)\ 

s is an indefinite fold point, or (iii) {u^y^ ^ ux — /?>)'. s is a 

cusp point, for some local coordinates (w, x, y) around s and an 
appropriate set of local coordinates around /(s) in the range M?. 
Moreover, no cusp point is a double point of / restricted to the 
set of singular points and / restricted to the set of all definite 
and indefinite fold points is an immersion with normal crossings. 
Thus for a proper stable map, a singular fiber-component passing 
through a definite fold point or a cusp point contains exactly one 
such point, while a singular fiber-component passing through an 
indefinite fold point may pass through one or two indefinite fold 
points. Otherwise, the map is called an unstable map. 

We note, characterizing the stability of the maps on compact 
3-manifold domains with boundary needs additional types of sin¬ 
gularities which is discussed in WT\ . Figure [Ta| is an example of 
a map from M? to where / = z). All its singular 

points are definite fold points, so this is an example of a stable 
map. Figure [Tc] is an example of a map from M? to M? where 
/ = -h+ z"^ — 5-1- z^) + 10, z). It has singular fiber- 
components which pass through four indefinite fold points (on 
corresponding four 1-manifold components numbered as 5, 6, 7 
and 8 in Figure O and so is an example of an unstable map. 

From the pre-image theorem (221, generically a regular 
fiber /“^(c) is a (J — r)-manifold for the regular value c = 
(ci,C 2 , ...,Cr). We note for d < r, f~^{c) is an empty set 
or a discrete set of points. A fiber f~^{c) can be consid¬ 
ered as the intersection of the fibers of the component scalar 
fields (ci), / 2 “^ (c 2 ),..., (cr) and a connected component 
of this intersection is a fiber-component. Alternatively, fiber- 
components of (/i, / 2 , ..., /r) can be considered as the contours 
of a component field f, restricted to the fiber-components of the 
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Figure 1: (a) A stable bivariate field (/i,/2) = y'^ , z) from M? to M? that is visualized using the transparent isosurfaces 

of the first component field; black curves are the fiber-components of the bivariate field; the red line represents the Jacobi Set; (d) 
The Reeb Space corresponding to (a) that is comprising one sheet (in pink) and the Jacobi Structure (red parabolic curve); (b) The 
Jacobi Set (consists of the red lines, top face and bottom face of the box) of the bivariate field (/i,/ 2 ) = y'^ z^, z) in the 

box [—1, 1] X [—1, 1] X [0, 1]; singular fibers passing through the boundary tangent points form a cylindrical surface that separates 
the domain into five components, denoted as A, B, C, D and E; (e) The Reeb Space of the multi-field corresponding to (b) that is 
comprising five sheets (in grey) and the Jacobi Structure (red lines); the regular components of the Reeb Space are marked to match 
the corresponding components in the domain; components of the Jacobi Set in the domain and their corresponding projections in 
the Reeb Space are denoted by numbers; (c) An unstable bivariate field (/i 5/2) = (v"^ — 5(x^ + z^) + 10, z) from M? 

to M? that is visualized using the transparent isosurfaces of the first component field; black curves are the fiber-components of the 
bivariate field; the Jacobi Set consists of 9 red lines; (f) The Reeb Space corresponding to (c) that is comprising six sheets and the 
Jacobi Structure (6 red lines). 


remaining component fields. This is a key observation, we use in 
building our Multi-Dimensional Reeb Graph data-structure. 

Jacobi Set. 

The compact J-manifold domain X (C M^) of the map / can be 
expressed as X = X° U 5X where X° denotes the interior (the set 
of interior points) of the domain and 5X denotes the boundary 
(the set of boundary points) of X. In case the domain X is without 
boundary, 5X = 0 and X = X°. 

Now the interior Jacobi Set of the map / : X ^ 
W is denoted by JJ^ and is defined by the set JJ^ := 


{x G X° I rank df° < min{J, r}} lIT^ where /° is the restriction 
of / to X°, i.e., /° := /|x° : -> W. In other words, JJ^ is the 

set of singular points of the map / interior to the domain X. Sim¬ 
ilarly, the boundary Jacobi Set of the map / is denoted by JJ^ and 
is defined as the set of singular points of the restriction of / to the 
boundary 5X, i.e. •= /|^x • Finally, by the Jacobi 

Set of the map / : X ^ we mean the union of the interior and 
the boundary Jacobi Set of the map /, and is denoted by JJy, i.e. 
J/ = J}UJ^. 

Now the boundary of a domain may come with corners, e.g. 
a 3-dimensional cube has corners of two types: 12 edge corners, 
and 8 vertex corners (as in Figure [T0). Let X be a compact 3- 
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dimensional manifold with corners and / : X ^ be a smooth 
map. A point q e dX is a. (boundary) regular point if / restricted 
to dX is a local homeomorphism around q, where dX stands for 
the boundary of X which includes all the boundary points and 
corner points. Otherwise, ^ is a (boundary) singular point. For 
example, if we take a point in a vertical edge in Figure [TH then 
in its ( 2 -dimensional) neighborhood on the boundary, there are 
always a pair of points that are mapped to the same point. Thus, 
it is never injective, and hence is never a local homeomorphism. 
Therefore, the point is a (boundary) singular point. 

Alternatively, the Jacobi Set is the set of critical points of one 
component field (say fi) of f restricted to the intersection of the 
level sets of the remaining component fields. Edelsbrunner et 
al. |[T^ studied properties of the Jacobi Set for r Morse func¬ 
tions. They proved the Jacobi Set is symmetric with respect to 
its component fields. They also showed, generically, the Jacobi 
Set of two Morse functions is a smoothly embedded 1-manifold 
where the gradients of the functions become parallel. However, 
in general Jacobi Sets are not sub-manifolds of the domain of the 
multi-field /, and are the disjoint union of sub-manifolds of the 
domain m. The red lines in Figures [Ta| and [Tel illustrate the 
Jacobi Sets of multi-fields on domains without boundary. 

Reeb Space. 

As with the Reeb Graph of a scalar field, the Reeb Space 
parametrizes the fiber-components of a multi-field and its topol¬ 
ogy is described by the standard quotient space topology. We 
note the fiber-components of a continuous map / : X ^ 

(X C M^) partition the domain X into a set of equivalence classes, 
denoted by W/ := X/ where two points a, Z? G X are equiv¬ 
alent or a b if f{a) = f{b) and a, b belong to the same fiber- 
component of f~^ and f~^ (/(^))- Now the canonical pro¬ 

jection map qf :X ^ X/ rsj that maps each element of X to its 
equivalence class defines the standard quotient topology where 
open sets are defined to be those sets of equivalence classes 
with an open pre-image, under map qf. The Reeb Space of / 
is the quotient space W/ together with this quotient topology. 
The decomposition of / as the composition of and /, where 
f \ W f is such that f = f oqf. This is called the Stein 

factorisation of /. The following commutative diagram describes 
this relationship between the maps. 

X- - -R'' 



Now to construct a fiber f~^ (a), instead of going directly from W 
to X one can compute the pre-image under f of a. Each fiber con¬ 
sists of a number of components, one for each point in 
Generically, the Reeb Space is a Hausdorff space, i.e., any two 
distinct points of W/ have disjoint neighbourhoods. Moreover, 
when r < J the Reeb Space corresponding to the multi-field / 
consists of a collection of r-manifolds glued together in compli¬ 
cated ways ca. 

Figures fTdl [Tel and fm show three examples of the Reeb Spaces 
corresponding to a stable bivariate field in M?, an unstable bi¬ 
variate field on a closed 3-dimensional interval and an unstable 


bivariate-field in M?, respectively. We indicate the dark (red) lines 
in the Reeb Spaces as the Jacobi Structures which are introduced 
in the next section. Note that the structures of the Reeb Spaces as 
in Figures [Tdl[Te| and [Tf| are obtained by analyzing the evolution 
of the fiber-components of the corresponding bivariate fields. For 
example, if we consider evolution of the fiber-components of the 
map in Figure [Tbl they start at the definite fold points on the line 
numbered as 1. Then these fiber-components start growing and 
meet at the boundary Jacobi Set points (on the lines numbered as 
2, 3, 4, 5). Then each of them splits into four fiber-components 
which continue to shrink and die at the corner singular points (on 
the lines numbered as 6 , 7, 8 , 9). This evolution phenomenon is 
captured in its Reeb Space (e). 

4 Theoretical Results 

In this section we exploit the underlying structure of the Reeb 
Space to decompose it into a set of simple manifold-like com¬ 
ponents (namely regular and singular components) and capture 
their connectivities by a dual skeleton graph (namely Reeb Skele¬ 
ton). Then, since in real applications most of the data come 
with simple domains (such as, a cube or a box), we study prop¬ 
erties of the Reeb Space and its representative skeleton graph 
for such topologically simple data-domains. More precisely, to 
prove our theoretical results in this section we consider a stable 
bivariate field / = (/i ,/2) : X ^ where X (C M^) is a three- 
dimensional bounded, closed interval. However, most of the re¬ 
sults are straight-forward to generalize for multi-fields of higher 
dimensions. Thus the domain X, we consider, is a compact do¬ 
main with boundary and is simply-connected and in this case the 
Reeb Space is path-connected. 

4.1 Path Connectedness 

For a continuous map / : X C ^ the Reeb Space is a quo¬ 
tient space of the fiber-components and is path-connected (or 0 - 
connected). That is, any two points po and pi of the Reeb Space 
can be connected by a path 7 : [0,1] -> W/ so that 7 ( 0 ) = po and 
7 ( 1 ) = pi. In other wards, we say 0-connectivity is preserved by 
the quotient map qf :X^Wf. This can also be stated by say¬ 
ing that 0-th homotopy group of the Reeb Space 7ro(W/) remains 
trivial. Next we prove the following important property of the 
Reeb Space. 

Lemma 4.1. Let f \ X C M? ^ be a continuous, generic 
map on a ?>-dimensional interval X and Wf be the corresponding 
Reeb Space. Let P be a continuous path between any two points 
on the Reeb Space. Then ifWf\Pispath-connected, then so is 

Proof. Consider any two points po^ p\ GWf\P. Since Wf\P 
is path-connected, 3 a path 7 : [0,1] -^Wf\P with 7 ( 0 ) = po 
and 7 ( 1 ) = pi. Using conditions like the genericity of /, Y{t) 
lifts to X \ qJ^{P), i.e., there exists a path 7 in X such that qyo 
f{t) = Y{t) (using the path-lifting property 1251 ). Now 7 is a 
path between any point of to any point of q'j^{p\) in 

X\q'j^{P). Therefore, X\qJ^(P) must be path-connected. □ 
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Thus Lemma l4Jl implies if there exists a path P in the Reeb Space 
whose preimage separates the domain then P must also 

separate the Reeb Space. This is a useful property in detaching 
unimportant components from the Reeb Space. 

4.2 Jacobi Structure 

As noted in Section[3l the Jacobi Set of a function is not the same 
as the set of singular fibers, as each point in the Jacobi Set is 
merely a representative of a singular fiber. Moreover, the struc¬ 
ture of the Reeb Space is actually given by a projection of the 
Jacobi Set or the singular fibers. For example, in Figures fTcl and 
m the Jacobi Set consists of 9 parallel lines in the domain, but 
they correspond to 6 1-manifold structures in the Reeb Space. 
Note that if the input multi-field domain is with boundary, there 
are additional edges (corresponding to the boundary Jacobi Set) 
needed to describe the Reeb Space. We therefore introduce the 
Jacobi Structure: the manifold structure of the Reeb Space cor¬ 
responding to the Jacobi Set in the domain. 

Definition 4.1. The Jacobi Structure of a Reeb Space W/ corre¬ 
sponding to a multi-field f : X CM? is denoted by Zf (^^d 

is defined by Zf := qf{Jf), i.e., the projection of the Jacobi Set 
Jif to the Reeb Space by the quotient map qf :X f. 

Note that according to our definition the Jacobi Set J/ consists of 
both the interior and the boundary Jacobi Set, i.e., J/ = JJ^ U J^. 
Thus each point of the Jacobi Structure corresponds to a singular 
fiber-component in the domain X of /, and vice-versa. 

To understand the underlying structure of the Reeb Space one 
needs to understand both the topology of the singular fibers and 
the corresponding local configurations of the Jacobi Structure. 
Classification of singular fibers and their local configurations in 
the quotient space have been studied for stable maps from M? to 
M? and to M? 1481 . Figure [2] illustrates examples of a regular 
and few singular fiber-components, and their local structures in 
the Reeb Space (341 [28l for a stable map / : X C ^ For 
a more complete classification of singular fibers for maps on 3- 
manifolds (with boundary) to plane and for local configurations 
of the Reeb spaces we refer to (4711461 . 

For a generic map / : X ^ W/ is a two-dimensional poly¬ 
hedron and Jacobi structure embedded in the Reeb Space consists 
of 1 -dimensional components which are at the boundary of the 
two-dimensional sheets in W/. Now a 1-manifold component 
of the Jacobi Structure can be classified into three types based 
on the transition of number of regular fiber-components if one 
passes across the component 113: 

1. Birth-Death or Boundary component - where a fiber- 
component takes birth or dies (Figured (b)), 

2. Merge-Split component or Bifurcation locus - where two (or 
more) fiber-components merge together or one component 
splits into two (or more) (Figure [2c)) and 

3. Neutral component - where there is no change in the number 
of fiber-components if one passes through such components 
(Figure Wig)), but here, the topology of the regular fiber- 
component changes from a circle to an arc (or vice versa). 



(a) regular point 


(b) definite fold 



(c) indefinite fold 



(e) double cone 


(g) boundary tangency point 


Figure 2: Regular and singular fiber-components and correspond¬ 
ing local configurations in the Reeb Space. The Jacobi Structures 
are in red lines in the Reeb Space (TTI . 


A connected component of the Jacobi structure may also consist 
of a composition of these three types, e.g. in Figure |2d) the 
Jacobi structure component consists of a boundary and a merge- 
split component connected at a discrete cusp point. In Figure[2e) 
four merge-split components are connected at a double point on 
the Jacobi Structure. 

Figures ns and [Tel respectively show an example of 8 1 - 
manifold components of the boundary Jacobi Set (red lines in 
the boundary of the domain) and their corresponding projection 
in the Reeb Space as 5 1-manifold parts of the Jacobi Structure. 
From this example, it is clear that a boundary Jacobi Set compo¬ 
nent may not be the boundary component of the Jacobi Structure 
in the Reeb Space or vice-versa. In Section [6l we propose an 
algorithm for computing the Jacobi Structure by constructing a 
Multi-Dimensional Reeb Graph corresponding to a multi-field. 


4.3 Regular and Singular Components 

As the number of dimensions increases, the projections of the 
singular fibers develop more internal structure in the Reeb Space. 
Consider the Reeb Graph of a scalar function: in this, the projec¬ 
tion images of the critical points are single points ( 0 -manifolds) 
separating edges (1-manifolds). Similarly, for the bivariate fields 
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shown in Figures [Tdl \T^ and [T2 the projections of the singu¬ 
lar fibers are arranged in a Reeb Space along 1-manifold curves 
which separate 2-manifold sheets. This induces a natural stratifi¬ 
cation or partition of the Reeb Space into disjoint subspaces (or 
strata). 

To describe a stratification of the Reeb Space and the corre¬ 
sponding domain of the multi-field we first classify the fiber- 
components of the generic map f :X CM? according to 

their complexity or codimension of the subspace where they lie 
Eg. Given the Stein factorization f = foqf, fiber-components 
of / can be classified into three classes. 

1. = {qj^{s) : 5’ G W/and does not contain any 

singular point of /}. Fiber-components of this class are the 
regular fiber-components and their ^y-images form codi¬ 
mension 0 subspaces in W/, denoted as W^. 

2. = {qj^{s) : s eWf and ^ (.s') contains exactly one de¬ 
finite or indefinite fold point}. Singular fiber-components 
of this class are moderately complex and their ^/-images 
form codimension 1 subspaces in W/, denoted as W^. 

3. = {qj^{s) : 5' G W/and ^}^^(5') contains a cusp point 
or two indefinite fold points}. Singular fiber-components of 
this class are the most complex and their ^y-images form 
codimension 2 subspaces in W/, denoted as W^. 

Complexity of a fiber-component increases as the codimension 
of the corresponding subspace in the Reeb Space increases. Note 
that ^y-images of the fiber-components in and form the Ja¬ 
cobi Structure 5/ of the Reeb Space, i.e., 3/ = Wj^ U W^. Topo¬ 
logically, regular fiber-components are either a circle or an arc 
ill. For stable maps / : X C ^ topologically there are 7 
different types of singular fibers in and 21 different types of 
singular fibers in Wf \ . 

Two regular points a^b E are topologically equivalent in 
the Reeb Space W/ or a Z? if there exists a path between a and 
h without intersecting the Jacobi Structure 3/. It is not difficult to 
check that ‘^p’ is an equivalence relation. Therefore, the equiv¬ 
alence relation ‘^p ’ partitions the regular points of W/ into a set 
of equivalence classes. Now we prove that each such equivalence 
class is a 2-dimensional sheet. 

Lemma 4.2 (Partition). The Jacobi structure d / of a Reeb space 
W f corresponding to a smooth stable map / : X C ^ sep¬ 
arates the Reeb Space into a set of 2-manifold components. 

Proof. Let Z) be a small disk in the range consisting of regular 
values (i.e., D does not intersect /(J/)). Then, by Ehresmann’s 
fibration theorem, / restricted to f~^ {D) is equivalent to the pro¬ 
jection D X F D, where E is a 1-dimensional compact man¬ 
ifold. So, this means that qf{f~^{D)) can be identified with a 
disjoint union of some copies of D, where the number of copies 
is the same as the number of connected components of F. Even 
when D intersects with /(Jy), if we restrict / to the components 
of the inverse image f~^{D) that do not intersect J/, then the 
same consequence holds. So, the regular sheets of W/ are lo¬ 
cally homeomorphic to D, and hence is a 2-manifold. □ 

Thus we have the following definition of regular components. 



Eigure 3: Reeb Skeletons: (a) corresponding to the Reeb Space 
in Eigure [Tdl (b) corresponding to the Reeb Space in Eigure [Tel 
(c) corresponding to the Reeb Space in Eigure [Tfl 

Definition 4.2. A path-connected component ofWf\Zf or 
is called a regular component. 

Generically, the 0-dimensional strata are in the boundary of the 
1-dimensional strata in the 5/- Therefore, an equivalence relation 
on the set of points in W| can be defined, similarly, where two 
points of Wj are equivalent if there exists a continuous path be¬ 
tween them without crossing the 0-dimensional strata in 5/ and 
each such equivalence class will be considered as a 1-singular 
component. 

Definition 4.3. A path-connected component of 23f\ or Wj- 

is called a 1 -singular component. 

Note that a 1-singular component in W/ may be an arc or a circle. 
An arc 1-singular component will also be called as an edge. 

Definition 4.4. Each component of is called a 0-singular 
component. 

To extract a skeleton graph from the Reeb Space we need ad¬ 
jacency of these regular and 1-singular components which are 
defined as follows. 

Definition 4.5. 1. A circle 1-singular component is self- 

adjacent (adjacent to itself). 

2. If two end points of an arc 1-singular component coincide, 
then the 1 -singular component is self-adjacent. 

3. Two distinct 1 -singular components 81^82 are adjacent if 3 
a 0-singular component Oq such that Si U ^2 U Ooform a 
connected space. 

Definition 4.6. A 1 -singular component Si is adjacent to a regu¬ 
lar component Rj if S ({JR j forms a connected space. 

Next we define a connectivity graph of regular and singular 
components based on their adjacency. 

4.4 Reeb Skeleton 

Once the Reeb space Wy is split into 2-manifold regular compo¬ 
nents and 1-or-lower manifold singular components, it is possible 
to perform a further reduction from the Reeb Space. To do so, we 
represent both these regular and 1-singular components as points 
(or nodes), and add edges representing their adjacency: in short. 
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(a) Reeb Space (b) Reeb Skeleton 

Figure 4: (a) Reeb Space with self-adjacent 1-singular compo¬ 
nent (b) Corresponding Reeb Skeleton. 

we can build the dual graph of these components of the Reeb 
Space. This has the merit of further reducing the Reeb Space 
from a 2 -dimensional structure to a fundamentally 1 -dimensional 
structure which is easier to represent, to reason about and to vi¬ 
sualise. We refer to this as the Re eh Skeleton and formally define 
as follows. 

Definition 4.7. Let R\^R 2 ^... ^Rmke the regular components and 
Si^S 2 ^ - • • ^Snbe the 1 -singular components ofW /. Then the Reeb 
Skeleton of f, denoted by Ky, is the adjacency graph which con¬ 
sists of (i) nodes nR. and nsj (i = 1 , 2 , ..., m and j = 1 , 2 ,...,^) 
corresponding to each of the regular and 1-singular components, 
and (ii) edges e{Sj,Sjf) and e{Ri,Sj) that are defined as follows: 

1. IfSj is self-adjacent, then e(Sj^Sj) = 1. In other words, nsj 
has a self-loop. 

2. If Sj is self-adjacent and Sj is adjacent with a regular com¬ 
ponent Ri, then e{Ri,Sj) = 2. In other words, nsj is con¬ 
nected with nR. by two edges. 

3. IfSj and Sjf are two distinct non-boundary I-singular com¬ 
ponents, then 

. x_ J 1 , if Sj^ Sf are adjacent 
I 0 , otherwise. 

4. For any regular component Ri and any I-singular compo¬ 
nent Sj 

( ^ \ _ j I 5 If^i^ Sj are adjacent 

I 0, otherwise. 

The regular and 1-singular components of the Reeb Space are 
represented as the regular and singular nodes, respectively, in the 
Reeb Skeleton. Figure [3] shows some examples of Reeb Skele¬ 
tons corresponding to the Reeb Spaces in Figure [T] Figure HI il¬ 
lustrates an example of the Reeb Skeleton with a self-adjacent 
singular node. Note that although the Reeb Skeleton gives a 
simple abstraction of 0-connectivity in the Reeb Space, it loses 
information of higher-dimensional connectivities, like higher di¬ 
mensional holes (tunnels, voids) in the Reeb Space. But on the 
other hand, the Reeb Skeleton is extremely useful for extracting 
any “fork”-like structure (corresponding to a merge-split feature) 
in the Reeb Space. And we will see later by a little simplifica¬ 
tion we can extract the most prominent merge-split feature in the 
Reeb Skeleton and so in the Reeb Space. Therefore, next we 
study properties of the Reeb Skeleton to simplify it further. 



Figure 5: Example of Reeb Space with a tunnel. 


4.5 Simple Domains 

We know from scalar fields that topologically simple domains 
have a useful property: the Reeb Graph is guaranteed to be a tree 
- i.e. the contour tree. This not only enables more efficient com¬ 
putation, but also provides straightforward mechanisms for fea¬ 
ture extraction, simplification and visualisation. Ideally, in multi¬ 
fields, the Reeb Space would also be contractible to a point. But 
we show this is not true, in general. 

In topology, simple domains are characterised by simply- 
connected space. A topological space is simply-connected if it 
is path-connected and every loop in that space can be continu¬ 
ously shrunk to a point without leaving the space. In terms of 
homotopy theory this means a simply-connected space is without 
any “handle-shaped hole” (as in Figure [5]) or it has trivial funda¬ 
mental group. For example, a sphere (that has a hollow center) 
is a simply-connected space whereas a torus (that has a handle¬ 
shaped hole) is not. Even a simpler topological space is known as 
contractible space which is homotopically equivalent to a point. 
Note that a contractible space is simply-connected, but the con¬ 
verse is not true. For example, a sphere is simply-connected as 
every loop on it can be contracted to a point on it, although the 
sphere is not a contractible space because of the center hole in 
it. In the following lemma, we prove that the Reeb Space cor¬ 
responding to a map defined on a simply-connected domain is 
simply-connected, but later we show it may not be contractible. 

Lemma 4.3 (Simply-Connected). The Reeb Space of a generic 
continuous map / : X C ^ A simply-connected. 

Proof. We consider any loop in the Reeb space W/. Then, it lifts 
to an arc in X. But, every fiber of is connected, and therefore, 
it lifts to a loop. As X is simply-connected, this lifted loop is null- 
homotopic. Therefore, its ^/-image is also null-homotopic from 
the continuity of This means that W/ is simply-connected. □ 

Therefore, if / is good enough (for example, triangulable or 
piecewise linear), then the Reeb space is simply-connected. This 
implies that the 1st homology of the Reeb Space also vanishes (or 
is the trivial group), and therefore the Reeb space does not have 
a tunnel or 1-dimensional hole (i.e., a hole inside a circle S^, e.g. 
Figure [ 5 ]). Thus we have the following theorem. 

Theorem 4.4. The Reeb Space of a generic map /: X C ^ 

does not contain any tunnel or 1-dimensional hole. 

On the other hand, for void or 2-dimensional hole (i.e., hole 
inside a sphere 5^), this is no longer true. We can construct 
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a (piecewise linear) map / : X ^ whose Reeb space does 
have a 2-dimensional hole. For example, consider the Hopf fi- 
bration ^ and its composition with a standard projection 
^ M?. The resulting map ^ M? is not generic, but per¬ 
turbing it slightly along its Jacobi set, we can obtain a generic 
map M?, whose Reeb space is the union of a 2-sphere and 

an annulus attached along the equator (and one boundary compo¬ 
nent of the annulus). Then, by extracting a 3-ball in the preimage 
of a two disk in the interior of the annulus part, we get the desired 
map X ^ M?. The Reeb space is the same space; the union of 

and an annulus (Figure©. Over each blue point lies a point 
(definite fold) and it corresponds to a birth-death. Over each red 
point lies a fiber as in Figure [2c) (with an indefinite fold) and 
the splitting of a circle fiber occurs. Over each green point lies a 
circle touching the boundary of the domain cube X. Thus, over 
each point in the shaded disk bounded by the green circle lies an 
interval. Note this disk is a subset of the annulus part. There¬ 
fore, a Reeb Space of a multi-field on a contractible domain may 
not be contractible and simplification of such space may not be 
simple as in the scalar case. 

According to Theorem l44l we can conclude that each regular 
component of W/ is planar; i.e., each regular component is a 
disk possibly with holes. For example, torus with holes (or a 1- 
dimensional hole as in Figure© never appears! This is essential 
in applying our simplification rules for the Reeb Skeleton as will 
be discussed in Section lb^l fFigurefTTI). Next we focus on finding 
a criterion for detachability of such regular components from the 
Reeb Space for simplifying the corresponding the multi-field. 

4.6 Detachability 

In the case of scalar field in a simply-connected domain, the Reeb 
Space (Graph) is a contour tree and there always exists a leaf edge 
that can be detached in a mathematically correct way, unless the 
contour tree consists only of one edge. We find similar criteria 
for defining detachable regular components in the Reeb Space. 

We say that it is possible to detach a regular component from 
a Reeb Space to obtain a simplified Reeb Space if the multi-field 
corresponding to the initial Reeb Space could be simplified to 
the multi-field corresponding to the modified one, and then the 
regular component is said to be detachable from the Reeb Space. 
Mathematically, any map could be simplified to a simpler map 
in the following sense. Since M? is contractible, any two stable 


maps /o and /i : X ^ are homotopic. So, using singularity 
theory, we can show that /o and /i are connected by a generic 
1-parameter family of maps. So, if we take an arbitrary stable 
map as /o and a very simple map as /i, then /o is simplified 
to fi after the generic 1-parameter family. Such a 1-parameter 
family passes through finitely many bifurcation parameters, and 
such bifurcations can be classified |[38l . 

Such transitions of the Reeb Spaces for generic smooth maps 
on a closed 3-dimensional manifold into M? have been studied in 
1^ . although for maps on a 3-dimensional manifold with bound¬ 
ary these results need further extension. In the current paper, we 
consider only a simple type of singularities and show that corre¬ 
sponding regular component is detachable from the Reeb Space. 
These components are known as lips and are defined as follows. 

Definition 4.8. A lip is a regular component that is attached to 
the other sheets of the Reeb Space exactly along one edge or an 
arc 1 -singular component, and it should not contain any vertex 
on the boundary, except for the two cuspidal points (Figure^a)). 

(a) (b) 

Figure 7: Lip simplification: (a) Reeb Space with a lip, (b) Sim¬ 
plified Reeb Space. 

Next we prove the following lemma to show that the underlying 
map corresponding to a lip can be simplified. 

Lemma 4.5. For a generic bivariate field / : X C ^ M?, a 
'dip ” can always be detached safely. 

Proof. If we have a lip in the Reeb Space, there are three pos¬ 
sibilities as shown in Figure [H That is, we may consider the 
map near the inverse image of the lip as a 1-parameter family of 
functions on a piece of surface. Let Sq be a “piece of surface” 
(a cylinder or a square as in Figure [8]), and f : M, r G I, be 

the 1-parameter family of height functions as in Figure© Then, 
the original map / is equivalent to the map {xf) {ft{x)f), 
X G ^o, ^ G I, around the inverse image of a neighborhood of the 
lip by The Figure ©presents the three such families of func¬ 

tions. As we can see easily, these can be eliminated continuously, 
by just shrinking the “time interval” for which a pair of critical 
points appear. □ 

Thus, we see that lips are detachable and they can be simplified 
as in Figure © Therefore, we get our simplification rule for 
detaching the lip components as follows. 

Simplification Rule: Let Rt be a detachable lip component of 
the Reeb Space W/. Then we simplify the Reeb Space by (i) delet¬ 
ing Ri with its adjacent boundary 1-singular component and (ii) 
converting the attached arc 1-singular component (merge-split) 
and two 0-singular components (cusp vertices) as regular. 

Next we discuss the Reeb Space (Skeleton) simplification 
based on the rule developed in this section. 
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□ D D D □ 

Figure 8: The behavior of the stable map near the pre-image of a 
neighborhood of the lip. The red lines indicate the boundary of 
the domain. 

5 Reeb Space Simplification and Mea¬ 
sures 

In the real multi-field data because of noise very often there are 
“lip”-like components which occlude the original feature cap¬ 
tured by the Reeb Space. Therefore it is important to simplify 
such components to understand the topology of the underlying 
data. Given that it is possible to detach such “lip”-like regular 
components from the Reeb Space, we follow a similar strategy to 
that used for the contour tree 0. There, a leaf edge was chosen 
for pruning and removed from the tree. If as a result a saddle 
point became regular (i.e. 1-manifold), it too was removed, sim¬ 
plifying the graph further. By tracking which leaves, saddles and 
edges are removed, the branch decomposition m then gave a 
natural simplification hierarchy for any ordering of leaves. 

In any Reeb space where Lemma [43] applies, we can use the 
same strategy, building a simplification hierarchy in the process. 
To do so, we simply choose a detachable component and remove 
it from the Reeb Space as described in the simplification rule of 
Section l431 We illustrate this process in Figure[T2l where we pro¬ 
gressively remove detachable regular components from the Reeb 
Space, reducing the Jacobi Structure accordingly as much as de¬ 
sired. As in leaf-pruning of contour trees, “lip”-simplification 
reduces the number of regular components in the Reeb Space by 
one each time, and also remove components of the Jacobi Struc¬ 
ture, guaranteeing that the number of steps required is linear in 
the number of regular components of the Reeb Space. Moreover, 
the editing operations to update the Reeb Space, Jacobi Structure 
and Reeb Skeleton are constant at every step, making the simpli¬ 
fication effectively linear (in the number of regular components) 
once the order of reduction is known. Therefore, we study differ¬ 
ent measures to associate with the regular components (nodes) of 
the Reeb Space (Skeleton). 

5.1 Range Measure 

In simplifying the contour tree, Reeb Graph and Morse-Smale 
Complex, simplification can be defined by cancelling pairs of 
critical points according to an ordering given by a filtration - i.e. 
a sequence by which simplices are added to a complex. For any 
given filtration, a unique ordering exists, and the persistence of 
a feature is defined by the distance in the filtration between the 


critical points defining the feature. 

For scalar data, however, the order in the filtration is dictated 
by the isovalues associated with each vertex of the simplex, with 
the result that persistence can also be formalised as the isovalue 
difference between the critical points that cancel each other. In 
multi-fields, the persistence of a feature gives rise to tuples rather 
than a single value O, which does not naturally give rise to a 
total ordering of the features. 

This is however, not the only way to define a simplification or¬ 
dering. Carr et al. 0 showed that pruning leaves individually 
could be ordered by geometric properties such as area, volume 
etc. of the features defined by the contour tree. In this model, 
persistence is the vertical height of a feature corresponding to 
a branch of the contour tree, and removing leaves can be done 
with simple queue-based processing. Recently, Duffy et al. ca 
demonstrated that many properties of isosurfaces in scalar and 
multi-fields relate to geometric measure theory. In this model, 
statistical and geometric properties of a function are measured by 
integration over the range. Following a similar approach we in¬ 
troduce a range measure for computing area of the regular com¬ 
ponents using the induced measure from the range to the Reeb 
Space. Note that, in general a regular component of a Reeb Space 
is projected to the range with multiplicities: i.e., this map is an 
immersion, but may not be injective. 

Consider for example the Reeb spaces shown in Figure [T] for 
bivariate volumetric maps. Mathematically, range measure of a 
regular component in the Reeb space W/ is defined as the area 
of the 2-dimensional sheets with respect to the measure induced 
from the usual area measure of the range Euclidean space. The 
range measure of each regular component in the Reeb space is a 
fixed scalar value. Thus, there is a unique induced ordering for 
simplification. If two components have identical range measure, 
some form of perturbation will be required to guarantee a strict 
ordering. 

5.2 Geometric Measures 

Similarly, it is also possible to compute geometric properties of 
the regular components, either in the domain, in the range, or 
in some combination of the two, using geometric measure the¬ 
ory. As with the contour tree 0 , obvious properties of interest 
include the measure of the region’s boundary in the domain (con¬ 
tour length in 2D, isosurface surface area in 3D), the measure of 
the region in the domain (area in 2D, volume in 3D), the measure 
of the function over the region (a generalisation of the volume in 
2D, hypervolume in 3D), and so forth. However, as in that work, 
rules will be needed in each case for combination of measure 
with parents in the simplification hierarchy based on the theory 
in Section 1431 

5.3 Summary of Theoretical Contributions 

We have now completed the theoretical groundwork for practical 
simplification algorithm of Reeb Spaces. In particular our theo¬ 
retical results could be summarised as follows. 

1. The Reeb Space consists of regular components correspond¬ 
ing to regions in the domain of the function, and singular 
components describing their relationships. 
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2. The Jacobi Set in the domain does not capture all of the 
structure of the singular components in the Reeb Space, and 
the Jacobi Structure is needed to do so. 

3. The Jacobi Structure of the Reeb Space can be used to fur¬ 
ther collapse the Reeb Space into the Reeb Skeleton. 

4. Multifields with topologically simple domains can be sim¬ 
plified using a variation on the leaf-pruning used for contour 
trees. 

5. A Reeb Space measure and other geometric measures are 
introduced to guide the Reeb Space simplification process. 

We now turn to the practical and algorithmic part of this paper: 
how to simplify the Joint Contour Net, an approximation of the 
Reeb Space. 

6 Algorithm: Simplifying the Joint Con¬ 
tour Net 

In this section, first we introduce the Joint Contour Net, a graph 
data-structure that approximates the Reeb Space. As described 
in Q, the Joint Contour Net is a quantized approximation of the 
Reeb Space. Therefore, to avoid having duplicate terminology 
we will use the same terminology for JCN as what we have de¬ 
veloped for the Reeb Space, namely, Jacobi Structure, Regular 
Component, Singular Components, Reeb Skeleton etc. 

6.1 Joint Contour Net 

The Joint Contour Net (JCN) 171 [TSl approximates the Reeb 
Space W/ of a multi-field / = (/i, /2, • • •, /r) 
in a J-dimensional interval X. Let / = (/i, /2, • • •, /r) : M ^ 
be a piecewise-linear (PL) approximation of / corresponding to a 
mesh M of X. The idea of computing the JCN is based on quanti¬ 
zation of the fiber-components of /. The JCN / with a quantiza¬ 
tion level (or level of resolution) thq is denoted as JCN{f ^mq), 
where mq refers to how fine the rectangular mesh for the range 
is. 

A quantized level set of fi at an isovalue h G Z/mg is de¬ 
noted by Qf^^{h) and is defined as: Qf^^{h) := {v G M : 
(^) round(mg^ (A)) =h}. A connected component of the quan¬ 
tized level set in the mesh is called a quantized contour or a con¬ 
tour slab. The part of the contour slab in a single cell of the mesh 
is called a contour fragment. 

Now the first step of the JCN algorithm constructs all the con¬ 
tour fragments corresponding to a quantization of each compo¬ 
nent field. In the second step, the joint contour fragments are 
computed by computing the intersections of these contour frag¬ 
ments for the component fields in a cell. The third step is to con¬ 
struct an adjacency graph of these joint contour fragments where 
a node in the graph corresponds to a joint contour fragment and 
there is an edge between two nodes if the corresponding joint 
contour fragments are adjacent. Finally, the JCN is obtained by 
collapsing the neighbouring redundant nodes with identical iso¬ 
values. Thus, each node in the JCN corresponds to a joint contour 
slab (or quantized fiber-component) and an edge represents the 
adjacency between two quantized fiber-components (with quan¬ 
tization level mo) of /. 




Figure 9: (top) The joint contour fragments and their ad¬ 
jacency graph for a PL-bivariate field defined by the values 
{(5,0), (0,0), (5,0), (3,2)} at the vertices of a mesh of two trian¬ 
gles. (middle) The Multi-Dimensional Reeb Graph constructed 
from the JCN. The critical nodes of the MDRG are the ‘red’ 
nodes which form the Jacobi Structure, (bottom) Corresponding 
Joint Contour Net, with critical nodes from the MDRG marked 
in colour. 

















Note that one can build a multi-resolution JCN by increasing 
or decreasing the quantization level using a scaling factor for the 
ranges of the component fields. An example of a small JCN is 
given in Figure [9j but we refer the interested reader to JTI for de¬ 
tails. The following lemma shows that in the limiting case, when 
the quantization level increases and the domain-mesh becomes 
more refined, then the JCN converges to the corresponding Reeb 
Space. 

Lemma 6.1 (Convergence). Let / : X c ^ d > r, be 

a tiangulable continuous multi-field with the Reeb space W/. 

Choose an increasing sequence of quantization levels {m^q} for 
f such that m^Q is an integer multiple of mg for each n and 
lim = oo. Furthermore, let be sequence of sufficiently 

fine meshes ofX such that is a refinement of M^-i for each n 
and lim d{Mn) = 0 , where d{Mn) stands for the maximum of the 

n^oo 

diameters of the cells ofMn. Finally, let f^^^ : Mn W be the PL 
map associated with f corresponding to the mesh Mn. Then the 
sequence {7CAr(/W , m^Q ) I converges to Wf. 

Proof. Hiratuka et al. (271 show for a PL map f \ A ^ B 
of a compact polyhedron A into another polyhedron B, if we 
subdivide the range polyhedron B appropriately, then A is sub¬ 
divided accordingly and the quotient map : A W/ to the 
Reeb Space W/ is triangulable with respect to the triangula¬ 
tions. In the proof, it is also shown that the inverse image by 
of a small regular neighborhood of a vertex v in W/ is al¬ 
ways a regular neighborhood of {qf)~^{v) in A. This implies 
that if the quantization level is high enough, then the quantized 
fiber-component is actually a regular neighborhood of the cen¬ 
tral fiber-component. Consequently, we have a natural embed¬ 
ding po : JCN^^ ^ Wy, where JCN^^ is the set of vertices of 

the Joint Contour Net JCN{f^^\m^Q) for sufficiently large n. 
(For each quantized fiber-component, associate the central fiber- 
component.) Furthermore, as is shown in (771 . this embedding 
preserves the adjacencies. This implies that the embedding po 
extends to an embedding p : JCN{f^^\m^Q) ^ Wy. Hence, the 

required result holds, since the triangulations of JCN{f^^\m^Q) 
and Wy becomes finer and finer as n increases. 

If / itself is not a PL map, then we can consider its triangula¬ 
tion g'.A^B and obtain the required result for the triangulation. 
As the domain X is compact, this implies the same consequence 
for the original map / as well. This completes the proof. □ 

Next we see that the simplification will have four stages: 1. 
extraction of the Jacobi Structure from the JCN, 2. computing 
regular and singular components for construction of the Reeb 
Skeleton, 3. computation of measures for each regular node in 
the Reeb Skeleton, and 4. simplification by pruning nodes corre¬ 
sponding to the “lip” components. In practice, the first stage is 
the most difficult - identifying the regular components, and this 
requires an intermediate data-structure, which we introduce now. 

6.2 Multi-Dimensional Reeb Graphs 

The first step in detecting and analysing the Jacobi Structure is to 
identify the nodes in the JCN that capture changes in the topology 


- i.e. the quantized representatives of the Jacobi Structure. To 
do so, we exploit a simple property of the JCN - that the slabs 
can be arranged hierarchically, with the levels of the hierarchy 
corresponding to the individual fields. At the highest level of the 
hierarchy, the slabs are only defined by field f \, and are therefore 
equivalent to interval volumes: as such, we can compute the Reeb 
Graph for field f\ (see Figure [9] (middle)). 


Algorithm 1 CreateReeb Graph (G,/;) 

Input: A subgraph G of JCN and a chosen field f 
Output: The Reeb Graph RG with respect to field f 
1 : Create Union-Find Structure UF for field f. 

2: For each adjacent gi,g2 € G with fi{gi) = Mgj), 

GFAddigugi) 

3 : for each component Q in UF do 
4 : Create a node nQ in RG 

5 : Map graph node-id(s) and field-values from G to nq 

6: end for 

7 : Order nodes , • • •, according to f field values. 

8 : for edge eie 2 in G do 

9 : if ^ 1,^2 ^ components Cj Ck and fi{ei ) 7 ^ fi{e2) then 

10 : Add edge e{ncj , ) in RG if not already present 

11: end if 

12: end for 
13 : return RG 


Each slab (i.e. interval volume) of /i can be broken up into 
smaller slabs with respect to field /2 in a similar way (which 
form a subgraph G in the JCN), and the Reeb graph for these 
slabs computed similarly, as shown in Algorithm [B Proceeding 
recursively, we then compute a hierarchy of Reeb graphs, each 
of which represents the internal topology of a slab of the parent 
Reeb Graph with respect to the child’s field. We call this hier¬ 
archy the Multi-Dimensional Reeb Graph or MDRG and denote 
this as My. 

Computing the MDRG is straightforward once the full JCN 
has been extracted: we start with the JCN and compute the Reeb 
Graph for property f\ by performing union-find processing over 
the nodes of the JCN. This breaks the JCN into subgraphs corre¬ 
sponding to slabs in the Reeb Graph of property f \. The MDRG 
for each subgraph is then computed recursively, and stored in the 
node of the parent Reeb Graph to which its slab corresponds. In 
the process, the slabs get separated out into smaller and smaller 
components. 

We state this as an algorithm in Algorithm [2] and illustrate with 
a bivariate field in Figure [9j This algorithm is stated recursively 
for simplicity, but can also be implemented with queue process¬ 
ing for speed. Moreover, the division of subgraphs at each level 
into slabs can be performed more efficiently by exploiting the 
connectivity already encoded in the JCN. 

The principal value of the MDRG is that every node of the JCN 
in the Jacobi Structure is guaranteed to be a critical node of the 
finest-resolution Reeb Graphs (denoted as the critical nodes of 
the MDRG). This immediately gives a method of computing the 
Jacobi Structure once the MDRG is known m. 
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Algorithm 2 MultiDimensionalReeb Graph (G,/;-,... ,/^) 
Input: Graph G, fields ffr Output: MDRG M 
1: if / < r then 

2 : Let R = CreateReebGraph(G, fi) 

3: Store R as root node of M 

4: for Each slab of R do 

5: Extract subgraph G^ of nodes of G belonging to in 

R 

6: Compute Ms = MultiDimensionalReebGraph(G^,//+i, 

••■Jr) 

7: Store Ms at node s of R 

8 : end for 

9 : return M 

10: else 

11: returnM = 0 

12: end if 


6.3 Jacobi Structure Extraction 

Since every node belonging to the Jacobi Structure is guaranteed 
to appear as a critical node of the lowest level of an MDRG, the 
initial stage in Jacobi Structure extraction is simply to mark these 
nodes. Unmarked nodes are then guaranteed to be regular, and 
can be collected into regular components. Once this has been 
done, any remaining nodes that are adjacent to each other and to 
the same set of regular components are identified, as these form 
a 1 -singular component between the regular components. 

The first stage of this can be seen in EigureO where the critical 
nodes of the lowest level of the MDRG together mark all of the 
Jacobi Structure nodes in the JCN (in colour). 

6.4 Reeb Skeleton Construction 

Once we have extracted the Jacobi Structure for the JCN, it is 
straightforward to compute the corresponding Reeb Skeleton by 
creating a single node for each regular or 1 -singular component, 
and connecting them using the adjacency of components in the 
Jacobi Structure. In Eigure 16.41 we see an example of the Reeb 
Skeleton construction for a volumetric bivariate field. We note, 
in this case, the Reeb Skeleton has no detachable lip-like regular 
node whereas the Reeb Skeleton Eigure [T 2 | (d) has such detach¬ 
able nodes. 

Now in the simplification algorithm, the order in which de¬ 
tachable Reeb Skeleton nodes are removed is determined by the 
metrics associated with those nodes. Computation of such met¬ 
rics are described next. 

6.5 Computing Simplification Metrics 

Our simplification algorithm can use any desired measure of im¬ 
portance for components of the Reeb space, including but not 
limited to: 

• Range measure. As described in Subsection 15.11 we can 
measure the size of the regular components by the induced 
measure of the range. This is easy to approximate - in this 
case, by the number of unique JCN slabs (i.e. pixels in the 
range) that map to a given regular component. 



Eigure 10: (a) Bivariate Eield: {x^ — z, +z^) in a 

box [—5, 5]^, the Ted’ components are the Jacobi Set, (b) the 
Joint Contour Net with the Jacobi Structure (in red), (c) Regular 
Components, (d) the Reeb Skeleton. 

• Surface area. A regular component of the Reeb space is 
separated from other regular components by one or more 
singular components in the Jacobi Structure. Since the regu¬ 
lar components correspond to features and the singular com¬ 
ponents to boundaries between features, we can associate 
the area of the bounding surface with the regular compo¬ 
nent for the purpose of simplification. Eor the JCN, we can 
approximate this with the surface area of the fragments ad¬ 
jacent to the bounding region. 

• Volume. Similarly, we can measure volume in the domain 
for each feature represented by a regular component, and ap¬ 
proximate it by summing the volume of fragments mapping 
to a given regular component. 

• Other measures. As shown by Duffy et al. Cl and 
Carr et al. 0 , almost any geometric or other property of 
features can be used for simplification provided that it is 
correctly approximated and suitable rules for composition 
during simplification are established. 

6.6 Simplifying the Reeb Skeleton 

The lip-simplification rule of the Reeb Space, as described in 
Section l4^ can be translated similarly in the corresponding Reeb 
Skeleton. We note, according to Theorem 4.4, each regular com¬ 
ponent of our Reeb Space is always a disk possibly with holes. 
Therefore, the degree 2 regular node B as in a Reeb Skeleton Eig- 
ure fTlT a) is always a lip-like detachable node. Eigure fTlT b) shows 
the simplified Reeb Skeleton after pruning the lip-like node and 
attached singular nodes. 

Einally, we give simplification strategies of the Reeb Skeleton, 
based on geometric and range measures of the components. The 
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Algorithm 3 SimplifyReebSpace 
I nput: JCN JCN 
Output: Reeb Skeleton Kf 

1: Build MDRG Mf and Jacobi Structure 5 f from JCN JCN . 

2: Partition JCN into disjoint regular components 
C = {Ri ,..., by deleting 5/ from JCN. 

3: Partition 5/ into disjoint 1 -singular components ,..., 5^} 

based on adjacency to regular components in C. 

4: Use adjacencies of {Ri,... ... ,5^} to construct K/ 

(following defintion l4.7b . 

5: Push “detachable” nodes of Kf on priority queue PQ with 
priority determined by geometric measures. 

6: while PQ not empty (or priority is below a threshold value) 

do 

7: Pop node r from queue 

8: Prune r from Ky 

9: end while 

10: return Simplified Reeb Skeleton Kf. 





(b) 

Figure 11: Simplification rule for detaching a lip-like node (reg¬ 
ular node B in (a)) from the Reeb skeleton. 


Figure 12: (Simplification Demo) (a) Original JCN/Reeb Space 
of bivariate field (Paraboloid, Height) (b) Jacobi Structure, (c) 
Regular components (d) Reeb Skeleton (‘blue’ corresponds to 
regular components and ‘red’ corresponds to adjacent 1-singular 
components) (e) Simplified JCN (f)-(g) Simplified Reeb Skeleton 
using range measure. 


Reeb skeleton simplification method simplifies the Reeb skeleton 
and the corresponding Reeb space given a threshold (between 
0 and 1) adapting the previous approaches in the literature 0 
[54ll . The threshold represents a “scale”, under which detachable 
regular nodes of the Reeb skeleton are considered as unimportant 
(noise). The threshold is expressed as a fraction of the range 
of the metric used. It can vary from 0 (no simplification) to 1 
(maximal simplification). 

Figure [T2| demonstrates the simplification of components from 
the Reeb Space, of an unstable bivariate volumetric data. We 
use range measure for ordering the components. We note, regu¬ 
lar nodes 4 and 2 in Figure [T2| are not strictly the lips according 
to our definition of lips. However, a perturbation can be applied 
first to convert such components to lips and then lip simplifica¬ 
tions can be applied. In Figure[T2|we apply our lip-simplification 
rule directly at the regular nodes 4 and 2, sequentially (node with 


smaller measure is pruned first). 

7 Implementation and Application 

We implement our Reeb Space or JCN simplification algorithm 
using the Visualization Toolkit (VTK) (U. Details of the JCN 
implementation can be found in ||7][T5l. Our Reeb Space sim¬ 
plification implementation takes the JCN (a vtkCraph structure) 
as input and builds four filters: (1) The first filter computes the 
Jacobi Structure by implementing the Multi-dimensional Reeb 
graph algorithm, (2) The second filter builds the Reeb Skeleton 
structure by partitioning the JCN, (3) The third filter implements 
persistence and geometric measures and (4) The fourth filter im¬ 
plements the simplification rules of the Reeb Skeleton. We use 
a vtkTree structure to store the Multi-Dimensional Reeb Craph 
(MDRG): Reeb graphs at each level of MDRG are stored in a 
vtkReebCraph structure. For capturing the Reeb skeleton we use 
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the vtkGraph structure. 


Figure [Tsl Note that the final simplification results a simple Y- 
fork where two regular nodes correspond to the separated nuclei, 
while the third represents the exterior. 


Table 2: Data Statistics 


datasets spatial-dimensions slab widths no. of nodes (JCN) no. of edges (JCN) 


(Circle, Line) 

(29, 29, 1) 

(1,1) 

500 

1057 

(Paraboloid, Height) 

(40,40, 40) 

(1,1) 

1260 

2383 

(Sphere, Height) 

(40,40, 40) 

(1,1) 

1308 

2428 

(Paraboloid, Sphere) 

(40,40, 40) 

(1,1) 

6554 

12795 

(Cubic, Height) 

(40,40, 40) 

(1,1) 

3149 

5928 


A force-directed graph-layout from the OGDF - an Open 
Graph Drawing Framework ca strategy has been used for the 
graph visualization as shown in the demonstrations and outputs. 
We run our implementation on different synthetic and simulated 
data sets for testing the performance. In Table[2]the synthetic data 
sets are labelled by the combination of scalar fields used: Circle: 

Line: y, Sphere: Paraboloid: — z, 

Height: z and Cubic: (y^ — xy + z^, v). Circle and Line are in the 
2D-box [— 5 , 5 ]^ and other fields are considered in the 3D-box 
[- 5 , 5 ] 3 . 


Table 3: Performance results for Simplification 


Data 

spatial 

Dimensions 

Slab 

Widths 

Jacobi 

Structure 

Reeb 

Skeleton 

Simplification 

(Circle, Line) 

(29, 29, 1) 

(1,1) 

0.06s 

0.242s 

0.00s 

(Paraboloid, Height) 

(40, 40, 40) 

(1,1) 

0.10s 

L97s 

0.00s 

(Paraboloid, Sphere) 

(40, 40, 40) 

(1,1) 

0.80s 

33.02s 

0.00s 

Nucleon 

(40, 40, 66) 

(8,2) 

0.45s 

48.91s 

0.45s 


Performance results Table [3] shows the performance results of 
the JCN and MDRG algorithms for some simulated data. All 
timings were performed on a 3.06 GHz 6-Core Intel Xeon with 
64GB memory, running CSX 10.8.5, and using VTK 5.10.1. 

The number of nodes in the MDRG is actually the number of 
Reeb graphs computed by the MDRG Algorithm [2l From the ta¬ 
ble it is clear that performance of the MDRG algorithm is quite 
impressive for these simulated data. The complexity of the Cre- 
ateReebGraph on a graph with n nodes is 0{n + plogn) which is 
the complexity of a sequence of p UF operations (here, p < n) 

lEl. 

Nuclear Scission Data 

In a previous application paper da, the JCN was applied to 
nuclear data set (time-varying bivariate field of proton and neu¬ 
tron densities) and used to visualise the scission points in high¬ 
dimensional parameter spaces. Here the scission refers to the 
point where a single plutonium nucleus breaks into two frag¬ 
ments. However, this was based on visual analysis, and was com¬ 
plicated by a number of artefacts such as the recurring chains of 
star-like motifs within the JCN. Moreover, the eight corners of 
the domain boundary induced eight corresponding small sets of 
features in the JCN. 

Here we apply our simplification algorithm to overcome these 
artefacts and preserve the principal topological feature. Note that 
apparently there are two red nodes that are adjacent to 5 blue 
nodes in the Reeb Skeleton Figure Ob). This never happens if 
the given multi-field is stable, so this is an example of an un¬ 
stable bivariate field in 3-dimensional interval. We demonstrate 
the process of simplification for one of these scission data sets in 


8 Conclusions and Discussions 

In this paper, we provide a rigorous mathematical and compu¬ 
tational foundation of multivariate simplification based on lip- 
pruning from its Reeb Space and this generalises approaches that 
are effective for scalar fields. We note, lip-simplification can be 
applied only when there is a lip-component in the Reeb Space 
and might not always be possible, but nevertheless, it is quite ef¬ 
fective when applied to real data sets, which usually contain a 
lof of noise (as demonstrated in Figure [13]). The Jacobi Structure 
that characterises the Reeb Space is richer than the Jacobi Set and 
decomposes the Reeb Space. This is proved to be a useful prop¬ 
erty for the simplification procedure. In addition, we have shown 
how to extract a reliable approximate Jacobi Structure and Reeb 
Skeleton from the JCN that can be simplified to improve the use 
of the JCN for multivariate analysis, and illustrated this with an¬ 
alytical datasets and a real-world data. However, there are few 
open issues which need to be addressed in future research. 

• False lips: Currently using our lip simplification approach 
we are not able to distinguish or simplify the false-lips, simi¬ 
lar as in FigureCll This is because our Reeb Skeleton cannot 
compute the multiplicity of the adjacency between a regular 
and a 1-singular node which will be important for detecting 
such a false lip in the Reeb Space. So, in our current simpli¬ 
fication we assume no false lip appears and we can simplify 
only the “genuine” lips as defined in Definition 14.81 



Figure 14: False lip 


• Discontinuity in components of Jacobi Structure: Com¬ 
puted 1-singular components of the Jacobi Structure in the 
JCN may be discontinuous because of the degeneracy, as all 
degenerate singular points may not be captured by the criti¬ 
cal nodes of MDRG. Moreover, a choice of the quantization 
level may also result in the discontinuous Jacobi Structure 
in the JCN. 

• Further Structures in the Reeb Skeleton: In the current 
implementation of the Reeb Skeleton we have considered 
only the regular and the 1-singular components and their 
adjacency graph. But, there is further hierarchy possible in 
the Jacobi Structure. For more than bivariate case, singu¬ 
lar components can be decomposed into lower dimensional 
manifolds (strata) and can be represented in the Reeb Skele¬ 
ton, hierarchically. However, detecting such lower dimen¬ 
sional strata needs further theoretical analysis and an algo¬ 
rithm for detecting them. 

Apart from these issues, in the future, we intend to work on fur¬ 
ther simplification and acceleration of these techniques, and on 
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Figure 13: This figure shows the simplification process of the JCN corresponding to a “nuclear scission” data set used in ifTSll . 
The process includes: (a) The original JCN graph; (b) The Reeb Skeleton; (c) A modified representation of the Reeb Skeleton by 
changing the colour of the degree 2 singular nodes as “blue”, since they capture only minor topological information, (d) Simplified 
Reeb Skeleton after applying lip-pruning using the range measure for ordering the nodes (e) Resultant Reeb Skeleton is represented 
as the “Y”-fork; (f) Geometry corresponding to the “scission” point. 


alternate methods for Reeb Space computation and / or approxi¬ 
mation. We also expect to examine more data sets from multiple 
domains, now that we have solved more of the main theoretical 
issues. 
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